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Abstract 



The meaning of 'tunneling' in a timeless theory such as quantum cosmology is dis- 
^Jjl cussed. A recent suggestion of 'tunneling' of the macroscopic universe at the classical 

turning point is analyzed in an anisotropic and inhomogeneous toy model. This 'inhomo- 
geneous tunneling' is a local process which cannot be interpreted as a tunneling of the 
universe. 



1. Timelessness and tunneling 

Quantum gravity and quantum cosmology have the reputation of not describing any accessi- 
ble experiments and for not having any consequences after the Planck era. However, it has 
been suggested that quantum cosmology could be responsible for some phenomena of everyday 
physics as the arrow of time and the classical appearance of our universe the smallness 

of the cosmological constant |[, or the onset of inflation p, |[. 

As a prominent subject over the years, tunneling processes have been considered for differ- 
ent purposes in quantum cosmology. In a recent paper D§browski and Larsen discussed the 
tunneling of a recollapsing universe at the classical turning point into an expanding state with 
increased scale factor Similar models have been studied e.g. also by and, for the early 
stage of the universe, by In this letter I will comment on the meaning of 'tunneling' in 

quantum cosmology and the model of |^ is substantially generalized. 

The Wheeler-DeWitt equation = which governs quantum cosmology may be regarded 
as analogous to the stationary Schrodinger equation. The similarity is particularly striking for 



the popular one dimensional models in which the universe is described by the scale factor 
a only. The universe then resembles a particle in one dimension and the Wheeler-DeWitt 
equation reads 



with a G 1R+. However, this similarity is misleading since, among other differences, there exists 
an external time parameter in ordinary quantum mechanics in contrast to quantum cosmology. 

The quantum mechanical time parameter is crucial even for the most simple examples 
of tunneling. Consider the reflection and tunneling of a wave function at a square potential 
barrier. This situation is described by superposing an ingoing wave with a reflected wave in one 
free region while in the free region on the other side of the potential barrier the wave function 
consist solely of an outgoing wave, the amplitude of which determines the tunneling probability. 
The concepts of in- or outgoing waves are justified by either considering a wave packet by 
superposing different energy eigenstates the sum of which resembles a particle moving in time. 
(Due to time-reparametrization invariance there is only the 'zero energy solution' to deal with 
in quantum cosmology. In one dimension it is thus impossible to form wave packets.) Or more 
simply by observing that the crests of the solitary waves are moving with ±/ca; — ujt = 2Txn. 
The necessity of a time parameter is, however, clear from the onset since the very notion of 
tunneling assumes a state changing in time: While initially there is no particle in some region 
(and classically there will never be one), there is a finite probability to find one subsequently. 

It is worthwhile to point out that a tunneling probability means probability of a tunneling 
'event' of a particle and thus implies a measurement. The concept of tunneling thus presupposes 
both a time parameter and a theory of measurement. Concerning the first point, it does not 
help that there is a conserved Klein-Gordon current in quantum cosmology in more than one 
dimension, since its sign cannot be fixed in the absence of an external time. Denoting something 
as outgoing as e.g. in the definition of the 'tunneling wave function' would completely 
arbitrary fix a direction on the configuration space. 

In the case of a high and broad potential barrier the situation simplifies, since the wave 
function in the forbidden region can be approximated by the exponentially suppressed solution 
only. In quantum mechanics the ratio of the (squared) wave function at the beginning and at 
the end of the forbidden region gives the tunneling probability. The same procedure is usually 
adopted in quantum cosmology as e.g. in in order to calculate the tunnehng probabihty of a 
bubble of false vacuum between two classically allowed regions in a 'free lunch' process. Accept- 
ing for the moment this interpretation, one arrives at similar conclusions in the semiclassical 
limit of quantum cosmology as in ordinary quantum mechanics. 

Usually, the situation is more complicated than in the examples above, as e.g. in the alpha 
decay when one of the two classically allowed regions is restricted to a finite region in space. 
But while in that case it is still possible to have a purely outgoing wave outside the nucleus, 
this breaks down if both classically allowed regions are finite. (D§browski and Larsen used this 
kind of potential for analyzing the tunneling probability at the classical turning point of a FRW 
universe with several matter sources.) In Fig. 1 the first situation is depicted by the dashed line 
while the case with two bounded regions is shown by the solid line. In order to calculate the 
tunneling probability it is thus not sufficient to consider stationary waves only. The calculation 
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Fig. 1: The solid line shows the kind of poten- 
tial used by D§,browski and Larsen. The sec- 
ond classically allowed region is due to domain 
walls and a negative cosmological constant. In 
order to get a 'tunneling universe', obviously, 
the ascending part of the potential is not nec- 
essary. In Sec. 2 a potential is used which is 
depicted by the dashed line. The decreasing 
is due to a positive cosmological constant. 



is much more involved as can be seen e.g. by the deviation from the exponential decay law for 
the nuclear decay in a box . Without the aid of the external time parameter the situation 
in quantum cosmology is quite unclear. One cannot even use the above mentioned comparison 
of the wave function at both ends of the potential barrier as a formal tool, since there will be 
no purely exponentially suppressed solutions. 

In view of the aforementioned problems I will use the notion 'tunneling' only as a formal 
notion provided there exists a purely exponentially decaying solution: The ratio of the squared 
wave function at the beginning and at the end of the exponential region is defined as 'tunneling 
probability'. (Note that the exponentially increasing wave function will also be a solution — 
'reversed tunneling' — and arbitrary superpositions of this two basic solutions.) 

This formal concept may only be regarded as a corresponding to a 'real process' if there 
is (at least) a time parameter and a theory of measurement in quantum cosmology. In more 
complicated examples as e.g. the one of Fig. 1 one furthermore has to be rather careful with 
the calculation of the probabilities. One may try to circumvent some of these difficulties 
by introducing a semiclassical time parameter e.g. due to a decoherence process or due to 
a Born-Oppenheimer type of approximation. Usually this point of view is put forward by 
mentioning that the tunneling occurs in a region where the semiclassical approximation should 
hold, although Kiefer and Zeh have argued that this argument is not sufficient |11|. If, however, 
decoherence could be used to define a semiclassical time, it is expected to simultaneously 
suppress the tunneling process by a Zeno-type effect ||12[. I will come back to these issues 
in Sec. 3 when the results of the calculations in the toy model are discussed. As frequently 
stressed above, the notion of tunneling makes sense only if there is a time parameter. Since 
this is certainly not the case in the Planck era one cannot sensibly speak of 'tunneling from 
nothing' of the universe which is frequently used as a quantum alternative to the big bang. 

The preceding considerations implicitly assume the so called naive interpretation of the 
wave function which itself relies on the similarity of the Wheeler-DeWitt equation with the 
Schrodinger equation. That is, I assume the wave functions to be normalizable on the whole 
(unconstrained) configuration space or that at least a conditional probability can be used. Any 
other interpretation in quantum cosmology which works on the reduced phase space or isolates 
a time parameter is meaningless in the example of the FRW universe with phenomenological 
matter. This is because in these interpretations there is but one physical state. But even for 
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more realistic examples in which the time parameter is given by some function of the volume, 
it does not make sense to compare volumes. As already mentioned above, even if a physical 
time has been introduced one has to establish a theory of measurement in quantum cosmology 
before one could sensibly speak of tunneling. 



2. An anisotropic and inhomogeneous model 

The Kantowski-Sachs model as the most simple anisotropic and even inhomogeneous model 
has the advantage that it can be solved exactly. The homogeneous version of the model com- 
bines spherical symmetry with a translational symmetry in the 'radial' direction. (First the 
homogeneous model is discusses before the translational symmetry is relaxed, see below.) The 
spacelike hypersurfaces of constant times are therefore cylinders 

els'" = dr^ + dQ^ . (2) 

Here b(t),z(t) G 1R+; b(t) is the surface measure of the two-spheres with metric dQ"^ and z(t) 
measures the spacelike distance between them; r is the radial coordinate. 

The model with positive cosmological constant A and pressureless dust is considered here. 
It is well known that it is only the presence of matter which renders the disklike singularity 
{z — > 0) from a mere coordinate singularity which indicates the incompleteness of the model into 
a curvature singularity. The dust is described by the parameter z,^ = pzb'^ = const (analogous to 
O'ui = PO'^ in the FRW model). For homogeneous models this approach is essentially equivalent to 
a more sophisticated one which starts from a Lagrangian for the dust degrees of freedom, see e.g. 
|[l3|l and the literature cited therein. Although dust as a matter source in quantum cosmology 
is unsatisfactory, it does here lead to a toy model which is calculable and complete. In addition, 
the present context of a macroscopic universe justifies this phenomeno logical description. 

The following account of the classical dynamics of this well known model is rather cursory. 
More details can be found e.g. in [|14| , |15| and the literature cited therein. That time parameter 



is used which is equivalent to conformal time in the FRW model (A^ = b). The dynamics is 
determined by the Hamiltonian constraint 



zb' + 2zbb + b^[z - z^- Azb') = , (3) 

(written in the configuration variables and their velocities) and by one of the equations of 
motion 

b^ + b-b,,,-^Ab' =: b^-VtJb) = 0. (4) 

Here the equation of motion for b{t) has already been integrated, with b^ as an arbitrary con- 
stant of motion. The equation of motion for the scale factor of the closed FRW model is identical 
to Eq. (^), d(t) = b(t), but the corresponding constant is fixed by the matter content. In 
the Kantowski-Sachs model the matter content instead fixes the constant of integration of z{t) 
as indicated by the notation. 

The classically forbidden regions in this model are determined by V;,m(^) ^ 0- Note that 
there is no forbidden region on configuration space due to the Hamiltonian (13) since its kinetic 
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term is indefinite. Consequently, the above defined forbidden region is defined for eacli classical 
solution separately, because each solution is uniquely represented by the 'effective mass' 6m- 
There are two classically allowed regions only if the cosmological constant is positive and smaller 
than the Einstein value: < A < 4/(96^). This potential which is depicted by the dashed line 
in Fig. 1 is similar to the situation in see the solid line in the same figure, except that it is 
not increasing for large values of b. 

The one-dimensionality of Eq. suggests, furthermore, the possibility of an inhomoge- 
neous z{t,r). Geometrically, this looks reasonable because a cylinder is defined only by the 
homogeneity of the surface measure b while an inhomogeneous spacelike distance z between 
the two-spheres would not deform it. It turns out that an inhomogeneous z{t,r) is even dy- 
namically consistent for pressureless dust as a matter content, as was first noticed by Ellis 
p6[] . A rigorous proof is possible by using the equations of motion for the general spherically 
symmetric model as is performed in the appendix. There it follows directly, since z' is always 
tied to b' and thus does not appear in the Kantowski-Sachs model. Due to this absence of any 
spatial derivatives, states at different points do not interact: The Hamiltonian (|^) remains thus 
essentially unchanged 

z{t,r)b^{t) + 2z(t,r)b(t)b(t) + b\t) {z{t,r) - z^{r) - Kz{t,r)b''{t)) =0, (5) 

One thus has a one parameter set of homogeneous solutions b = b{t) which then determines 
uniquely the solution z = z{t, r) the inhomogeneity of which is due to the inhomogeneity of the 
dust. 

In the recoUapsing region one gets a qualitatively correct picture of the dynamics by con- 
sidering the exact solutions for A = 0: 

bit) = b^sm^ and z{t,r) = K{r) cot (0 + ^^^(r) [l - ^ cot Q)) , (6) 

with t G [0, 27r] and where i^' is a physically meaningless 'constant' of time-integration (different 
functions of K are identified by a redefinition of the radial coordinate and by simultaneously 
changing the dust potential). The other classically allowed region for large b can be approxi- 
mated by considering only the cosmological term in the potential. Obviously, the inhomogeneity 
of z{t, r) is due to inhomogeneous dust z^{r). 

As usual Dirac's quantization scheme is used by turning the variables 6, z and their conjugate 
momenta into operators which satisfy the standard commutation relations. The Hamiltonian 
constraint (|^) is turned into an operator which annihilates the physical states: = 0, the so 
called Wheeler-DeWitt equation. In the configuration space representation it reads 



- 2zb^H'^ 



z^d^ + kzd, - 2zbdbd, + z^b^ (l - Ab^) - zz^,b^] 2) = , (7) 



where and are the partial derivatives with respect to z and b, respectively. Factor ordering 
is partially left open as indicated by the parameter k. The Laplace-Beltrami ordering is given 
by = 1. This equation can exactly be solved for a more general potential which contains 
arbitrary functions of b multiplied by z and 2;^ ^4 . 
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Although z{r) is a field it suffices to consider the minisuperspace Wheeler-DeWitt equation 
at each point r since there are no partial derivatives with respect to r (and because of the 
simple structure of the solutions). Quantum ffuctuation are expected to result in interaction 
between different points. But this effect can only be considered in the context of a more general 
model as is common to this kind of problems. 

In order to solve the Wheeler-DeWitt equation it is convenient to introduce the operator 
6in := —\dz^ + {h — Mr" /?)) since [6m, H] = 0. As indicated by the notation, the eigenvalues of 
this operator is the effective mass for b. The eigenvalue equation 6m^f)m = ^m^bm can easily be 
solved, and one finally ends up with the following set of exact solutions for the Wheeler-DeWitt 



equation 



\ 



db' 



(8) 



Alternatively, one can consider the solutions which are given by a superposition of the expo- 
nential with plus and minus sign to form e.g. the 'cos' and the 'sin' (the 'cosh' and 'sinh' in 
the forbidden region). The integral in the exponential can be expressed in terms of elementary 
function for 6m = (that is for large values of b) and for vanishing A (that is for small values 
of b). Neither case is appropriate here. 

In order to analyze the formal tunnel probability one has to calculate the ratio prob : = 
l^(^i)P/l^(^2)P where bi and 62 denote the beginning and the end of the classically forbidden 
region, respectively. While, the prefactor of the wave function is divergent at the borderline 
of the classically allowed region it cancels in prob, since Vb^{b) possesses three distinct, simple 
zeros (the third one is negative and has no physical significance). The first term in the exponent 
vanishes for b = 61,62- There remains the second term in the exponential which is a definite 
integral between 61 and 62. This integral remains finite since at both limits of integration the 
integrand is approximately 1/x and the infinities at the two boundaries cancel each other. 

If one is interested in a 'tunneling' from the recoUapsing region into the region right 
from the potential barrier, one has to choose the minus sign in the exponent of the wave 



function. (Note that the Hartle-Hawking wave function has the opposite sign |T5[ and de- 
scribes thus a 'tunneling' into the recoUapsing region.) The final result then is given by: 
prob = (61/62)^'^"^'' exp{— 2;m-F(6m, A)}, where the number F(6m,A) is the definite integral. In- 
serting some simple values: 6m = l,A = 3/8 gives 61 = y/E— 1 and 62 = 2; the evaluation of the 
integral gives: F(6m, A) pa 5.106. Choosing furthermore the Laplace-Beltrami factor ordering, 
k = 1, one finds 

prob ^ exp{— 5.1062;m(?")} • (9) 

The simple structure of this result relies on the semiclassical structure of the exact wave function 
plus the canceling of the prefactor in prob. 



3. Results and remarks 

A remarkable feature of the above 'tunneling probability' is the way it depends on the matter 
content: The logarithm of prob depends linearly on Zm- Apart from factor ordering effects the 



6 



tunneling probability equals one for the vacuum modelQ. Since the matter does not influence the 
tunneling barrier one has the analogue of a particle with mass (squared) running against 
a potential barrier. However, different points r in the Kantowski-Sachs model will 'tunnel' 
individually since they do not interact. This situation is analogues to a cloud of non-interacting 
particles. In more realistic models one might think of weakly coupled points (perhaps galaxies), 
which nevertheless will behave mainly independent in the 'tunneling process'. This local process 
clearly cannot be interpreted as a tunneling of a whole universe. 

How is this process to be interpreted? Due to the timeless nature of quantum cosmology, 
and in particular since a semiclassical time might not be defined at the turning point, both 
'universes' (that is both classically allowed regions) exist 'simultaneously'. In this case one might 
interpret the tunneling as a quantum wormhole^. However, the same argument of timelessness 
tells one, that there is no classical observer at the turning point, no incoming and no outgoing 
wave. The world is completely quantum and consequently one cannot speak of a tunneling 
process or of wormholes. If, on the contrary, the tunneling occurs as a change 'in time', 
e.g. defined by decoherence, the variable b in some space regions changes from femax of the 
recoUapsing solution to bmm of the solution right from the potential barrier, with fomin > &max- 
This is similar to the 'free lunch' process in the very early universe p, ^. These tunneling 
regions behave dynamically different from their environment since they are now expanding 
instead of recoUapsing and because of the sudden change of the 6-variable. This results in a 
destruction of the cylinder geometry of the universe. 

It has been emphasized by Dabrowski and Larsen that this tunneling process does not 
require any change in the matter content as e.g. a change of the vacuum of the involved fields. 
However, Rubakov has suggested that the matter content may be changed due to the very 
tunnehng process He considered a scalar field, conformally coupled to the FRW model, and 
observed that the tunneling probability is enhanced with growing particle content. Probably, 
this remark will apply too when only small parts of the universe are involved. However, the 
calculation for the Kantowski-Sachs model is much more involved than that for the FRW model. 

A comparison of prob for the Kantowski-Sachs and for the FRW model shows another 
difference than that of inhomogeneity. Starting from the Wheeler-DeWitt equation for the 
FRW model [dl + |Aa^ + a'^ — 0^0^]'^ = one gets (in WKB approximation) 

P™^FRw = exp |-2y" a^V{a,a^,A)da^ , (10) 

where V is defined as in the Kantowski-Sachs model. This result is quite different from the 
expression (^ for the Kantowski-Sachs model, since in the FRW model Om is determined by 
the matter content. The analogue for the tunneling in the FRW model is that of a particle 
running against a barrier the width of which is fixed by the matter content. 



It was shown in WM that one can insert a Kantowski-Sachs cylinder between two FRW 



half-spheres. The 'tunneling probability' in this compact model changes drastically at the 



^ In this way the model automatically takes into account that the vacuum model does not possess forbidden 
regions: It describes the dynamical regions of a black hole (in de Sitter space). 

^ A quantum wormhole is formally defined by an Euclidean region between two separated classical regions. 
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borderline between the different parts. However, since the dust content in the FRW model 
fixes the constant of integration fom of the Kantowski-Sachs region, the tunneling probability 
in the Kantowski-Sachs region is a function of the FRW dust, too. The functional dependence 
is even similar in both regions. 

It has been argued in the last section that the calculation of prob is not spoiled by the 
divergences of the wave function (P). One can even get rid of these divergences by considering 
superpositions. (In contrast to the delta-functional like solutions (H) it is e.g. possible to 



consider 'free waves', the form of which reveals nothing of classically forbidden regions |17].) 
This is to be done in any case since usually a universe is not represented by one eigenfunction 
but by a wave packet. The above calculation of prob is nonetheless necessary since the forbidden 
regions are defined for fixed values of b^ only. With other words, every single component of 
the superposition 'tunnels' independently. Moreover, in more complicated models there might 
not be a sharply peaked wave packet at the classical turning point due to interference effects 



between the 'incoming' and the 'outgoing' part |T^. This is in agreement with the point of 
view that there is no tunneling because the wave function at this point is completely quantum. 

There might be thus two effects of quantum cosmology at the classical turning point of 
the universe: The breakdown of the semiclassical approximation and a local change into an 
expanding state. One might try to circumvent the breakdown of the semiclassical approximation 
by considering decoherence. If this worked (Kiefer and Zeh have expressed their doubts that 
it does [0), the increased classicallity of the universe would further suppress the tunneling 
probability |jl2|. Both quantum effects would then be lost simultaneously, by the same token. 



Acknowledgment 

I thank Dieter Zeh and Glaus Kiefer for their enlightening remarks, Hans Kastrup for his 
criticism, and Henning Wissowski and Thomas Strobl for supporting me in the difficult task of 
making this text intelligible. This work was supported by the DFG. 



Appendix 

In order to proof that an inhomogeneous z{t, r) is possible in the Kantowski-Sachs model with 
dust one has to discuss the equations of motion of the general spherical symmetric model. Its 
metric can be parametrized by 

ds^ = -N\t, r)df + L^{t, r)dr^ + R^{t, r)dn^ (11) 

where a vanishing shift-function has been chosen. One nevertheless has to satisfy the super- 
momentum constraint 

Here and in the following all equations are presented with the fields and their velocities instead 
of their momenta. According to this constraint a homogeneous surface measure R{t, r) — »• 
R{t) = b{t) necessitates a homogeneous lapse function N{t,r) N{t). In the following the 
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gauge fixing N — R similar to the one in the main text is chosen, but the homogeneity of N, R 
is not yet enforced. 

The gravitational Hamiltonian reads 

° ^ + + ^^^^ 

and the equation of motion are given by 

" R^ R 3RR'^ A , .s 

[LRf = i(M) -^i! + fl^(|y + ALiiVfl(:^y . (15) 

The homogeneity requirement R(t, r) — > R{t) = b(t) obviously leads to an equation of motion 
for b{t) which is independent of any other variable. One can thus insert a solution of it into one 
of the other equations in order to determine L{t,r). Furthermore, since L' does only appear 
in combination with R' all partial derivatives cancel and the equations reduce to that of the 
inhomogeneous Kantowski-Sachs spacetime. 

A non-vanishing pressure of a matter source would lead to an interaction of neighboring 
points. Thus, only incoherent or pressureless dust might serve as a matter source in this 
model. Since according to the Bianchi identities the dust flow is geodesic one may consider 
radial moving dust. Furthermore, by the other part of the Bianchi identities one obtains matter 
conservation p{t, r)L{t, r)R'^{t, r) = Cm(r). This general case is known as the Tolman model and 
in the case R{t, r) — > b{t) one gets the inhomogeneous Kantowski-Sachs model. In the above 
equations there will be an additional term in the Hamiltonian (the potential is supplemented 
by — Cm(r)) but the equations of motion remain unchanged by the dust. Obviously, there are 
no further complications to the inhomogeneity argument due to dust. 
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